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The isospin diffusion and other irreversible phenomena are discussed for a two- 
component nuclear Fermi system. The set of Boltzmann transport equations, such 
as employed for reactions, are linearized, for weak deviations of a system from unifor- 
mity, in order to arrive at nonreversible fluxes linear in the nonuniformities. Besides 
the diffusion driven by a concentration gradient, also the diffusion driven by temper- 
ature and pressure gradients is considered. Diffusivity, conductivity, heat conduction 
and shear viscosity coefficients are formally expressed in terms of the responses of 
distribution functions to the nonuniformities. The linearized Boltzmann-equation 
set is solved, under the approximation of constant form-factors in the distribution- 
function responses, to find concrete expressions for the transport coefficients in terms 
of weighted collision integrals. The coefficients are calculated numerically for nuclear 
matter, using experimental nucleon-nucleon cross sections. The isospin diffusivity 
is inversely proportional to the neutron-proton cross section and is also sensitive to 
the symmetry energy. At low temperatures in symmetric matter, the diffusivity is 
directly proportional to the symmetry energy. 

PACS numbers: 21.65.+f,25.70.-z,25.75.-q 



I. INTRODUCTION 

The availabihty of beams largely differing in isospin content in nuclear reactions has 
dramatically increased interest in phenomena associated with the variation of that content. 
In the context of peripheral reactions, this includes interest in changes of the nuclear den- 
sity profiles with the isospin content. In central reactions, the attention has focussed, in 
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particular, on the dependence of isospin symmetry energy on density. The determination 
of that dependence would permit an extrapolation of the nuclear equation of state to the 
neutron matter limit The isospin asymmetry has been used in central reactions for 
projectile-target tagging in the investigation of stopping P]. 

This paper deals with the irreversible transport of isospin and other quantities in a 
nuclear system, as pertinent for reactions, for small deviations from equilibrium. In that 
limit, the irreversible transport acquires universal features and is characterized in terms 
of transport coefficients, that include the isospin-diffusion coefficients. The coefficients are 



derived here for the c 
reaction simulations 



namics described in terms of a Boltzmann equation set such as used in 
1, 0] . The main diffusion coefficient or diffusivity, characterizing isospin 
diffusion driven by the gradient of asymmetry, is evaluated using free neutron-proton cross 
sections. In the past, other transport coefficients, viscosity and heat conductivity, have been 
investigated for nuclear matter 0, Q, ^ O]- It was subsequently found that conclusions 
from comparisons of reaction simulations to data on stopping can be universally formulated 
in terms of the nuclear viscosity Q| . It is hoped that the diffusivity can be of such utility 
as that other coefficient, for the systems with a varying isospin content. 

The past studies of irreversible linear transport for nuclear matter were primarily directed 
at momentum and energy. Tomonaga '5l] and Galitskii et al. obtained the low- and high- 
temperature limits for the shear viscosity and heat conductivity. Danielewicz derived 
results for those coefficients valid in a wide range of nuclear densities and temperatures. 
Hakim and Mornas Q studied different transport coefficients within the Walecka model 
following the relaxation-time approximation. 

Our derivation of diffusion coefficients follows the general strategy of Chapman and En- 
skog , but here for a Fermi system, with inclusion of mean-field effects such as appropriate 
for a nuclear system. In the next section, we discuss the diffusion coefficient concept qualita- 
tively and make simple estimates for the nuclear matter. The modification of the Boltzmann 
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equation to extend it to fermions has been first discussed by Uhlenbeck and Uehling 
In Sec. Ill, we formally solve the set of Boltzmann equations for a binary system of fermions 
to find thermodynamic fluxes driven by speciflc thermodynamic forces and to flnd general 
but formal expressions for the diffusion and other transport coefficients. The transport coef- 
ficients have been (as we found) flrst considered for fermions by Hellund and Uhlenbeck ^ | ; 
compared to them, our notation here adheres more to what is now customary for nuclear 
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reactions. Closely related to the diffusivity is the electrical conductivity that is included 
in our considerations. In Sec. IV, we obtain more specific results for the coefficients on 
assuming deviations from equilibrium suggested by the Boltzmann equation set, for specific 
thermodynamic forces present. Numerical results for the coefficients are obtained in Sec. 
V, using free NN cross sections. We also estimate there the pace of isospin equilibration in 
reactions. We summarize our results in Sec. VI. More technical mathematical details and 
some reference information are provided in five appendices. In sequence, these appendices 
are devoted to the definitions of macroscopic quantities, the continuity equations, the conti- 
nuity equations for an ideal fiuid, the transformations in the driving force for diffusion and 
to the algebra of collision brackets. 

II. DIFFUSION IN A BINARY SYSTEM 

Diffusion and other irreversible transport processes occur when a system is brought out of 
equilibrium. The direction of those processes is to bring the system back to the equilibrium. 
For small perturbations, in terms of constraints that may be set externally, the system 
response is linear in the perturbation. The coefficient of proportionality between the induced 
flux and the perturbation is the transport coefficient. 

In a multicomponent system with no net mass flow, irreversible particle flows result if 
particle concentrations are nonuniform. For N components, there are N — 1 independent 
flows and N — 1 independent concentrations (since the concentrations need to sum up to 1). 
The flows are then related to the gradients of the concentrations with an {N — 1) x {N — 1) 
matrix of diffusion coefficients. In a binary system, only a single coefficient of diffusion, 
or diffusivity, relates the irreversible particle flow to the nonuniformity in concentration. 
However, as we shall see, nonuniformities in other quantities than concentration, can induce 
a dissipative particle ffow as well. 

Our focus, obviously, is the binary system of neutrons and protons. However, for the 
sake of utility of the results elsewhere and for the ability to examine various limits, we 
shall consider a general two-component system of fermions. An extension of those results to 
bosons, outside of a condensation, will be trivial. 

The two components will be denoted 1 and 2. Then, for the particle i, the density 
is rii — Ni/V, where Ni is the particle number in some inflnitesimal volume V. With 
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net density n — rii + n2, the particle concentration for 1 is = ni/n and for 2 it is 
l — u — U'i/n. Moreover, with representing the mass of particle i, the net mass density is 
p = P1+P2 = TTii ni+m2 n2, and the mass concentration for i is Cj = rrii rii/p. The differential 
particle concentration is 5 = {rii — n2)/n. The different concentrations are obviously related 
and thus we have v = {1 + 5)12 and ci = mi (1 + 5)/ (mi {1 + 5) + 1712 (1 — S)). Later in the 
paper, we shall primarily use the differential concentration S as an independent variable. 
The dissipative particle flows Fj are defined relative to the local mass velocity v, 

= n,(Y,-v), (1) 

where Vj is the local velocity of i'th component and 

V = (piXi + p2V2)/p. (2) 

We might consider other flows such as defined relative to the local particle velocity, but 
those flows are combinations of Fi and F2. Moreover, even Fi and F2 are redundant and 
we might just use Fi as an independent flow with the flow of 2, as easily seen, given by 
F2 = — miFi/m2. Another option might be to use as independent the differential flow 
defined as 

F^ = Fi - F2 . (3) 

If the system is at uniform pressure and temperature, but there is a small concentration 
gradient present, the fluxes develop linear in the gradient, enabling us to write, e.g. 

r. = -„D* aud r,^-nD,?^. (4) 
or or 

These are so-called Pick's laws. Notably, the stability of an equihbrium state requires > 0. 
Since mi Fi + m2 F2 = 0, we have mi Di = 1712 D2. For the differential flow, we have 

du (9(1 — u) ^ dS 

Ts = -nD^— + nD2 ^ ^ ^ ^-nDs — . 5 
or or or 

Here, the differential coefficient is Ds — {Di + D2)/2. 

So far, we assumed a system at a uniform pressure and temperature, with just con- 
centration changing with position. If the variations in a system are more complex, other 
nonequilibrium forces than the concentration gradient can drive the diffusion. This will be 
explored later in the paper. General guidance regarding the forces which can contribute 
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is provided by the Curie principle. This principle exploits symmetry and states that the 
driving forces must have the same tensor rank and parity as the flux they generate. 

For the system of neutrons and protons, the differential concentration 6 becomes a concen- 
tration of the isospin and the differential flow becomes the isospin flow, = F/. Moreover, 
the differential diffusion coefficient becomes an isospin diffusion coefficient, Ds = Dj, and 
for equal masses we expect Dj = Dp = D^. 

It is popular to relate the concept of a diffusion coefficient to a diffusion equation. Indeed, 
if we consider a uniform system of protons and neutrons at rest, but with the nucleon 
concentration changing in space, then, from the continuity equation for the differential 
density 

(^) 

we get the familiar equation 

f,-D.VH. (7) 

Here, for D/, we have assumed a weak dependence on the concentration 5. 

Before turning to a derivation of rigorous results for the diffusion and other transport 
coefficients, it may be instructive to produce simple mean-free-path estimates for those 
coefficients. Let us consider components of equal mass (the mass then becomes a simple 
normalization coefficient in density that may be factored out) and consider the gradient of 
concentration along the x axis, in the medium at rest. If we take the three coordinate axes, 
then 1/6 of all particles will be primarily moving along one of those axes in the positive 



or negative direction, with an average thermal velocity V_ = y^3T/m, for the distance of 
the order of one mean free path A, without a collision. Considering the particles 1 moving 
through the plane at x = 0, they will be reflecting density at a distance A away. Including 
the particles moving up and down through the plane, we find for the flux 

Ti^l {n,{x - A) - n,{x + \)) 'l^^^ ' 
With Q, we then get for the diffusion coefficient 



1 I T 

Dr^-XV J—, (8) 

d na \ om 

with A ~ l/(r7,o"). A more thorough investigation shows that it is the cross section for 
interaction between the two species that enters the diffusion coefficient. 
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Let us now evaluate the magnitude of the isospin diffusion coefficient. At temperature 



T ~ 60 MeV and normal density uq = 0.16 fm~^, with a^p ~ 40 mb, we find Dj ~ 0.2 fmc. 
We will see this to be in a rough agreement with thorough calculations. 

Similarly to the above, one could employ the mean-free path arguments to determine 
the better investigated coefficients: shear-viscosity rj and heat conduction k. One finds 
T] ~ ^nmV_X and k ~ ^nV_Xcv, where cy is the specific heat per particle. For T ~ 60 
MeV and a ~ 40 mb, we find r/ ~ 30 MeV/(fm^ c) and n ~ 0.06 c/fm^. Up to factors, the 
shear viscosity and heat conduction coefficients play the role of diffusion coefficients in the 
diffusion equation for velocity vorticity and in the heat conduction Fourier equation identical 
in form to the diffusion equation. 

In the estimates above, we just considered the free motion of particles in-between colli- 
sions. If self-consistent mean fields produced by the particles depend on concentration, then 
this dependence, on its own, contributes to the diffusion. In the case of nuclear matter, the 
interaction energy per nucleon may be well approximated in the form quadratic in isospin 
asymmetry, Ej = a'"j6'^, where 6 = {up — nn)/n and is the interaction contribution to 
the asymmetry coefficient ai. At normal density, the coefficient is ~ 14 MeV. The naive 
expectation for two-body interactions is that is linear in density. At constant net density, 
the quadratic dependence of the interaction energy on 6 leads to the force Fp^„ = =F-^ 
of opposite sign on protons and neutrons. The direction of the force for positive is to 
reduce nonuniformity in isospin. Under the influence of this force, a proton accelerates for 
a typical time between collisions At = \/V_ and then, in a collision, resets its velocity. The 
described polarization effect augments then the proton flow by 

Ar, = „,Av, = -l^-^^. (9) 

^ n 2my_ or 

In comparing with (jHl), after correcting for the local center of mass motion, we find that the 
polarization increases the diffusion coefficient by 

d;~(i-5^)|d?, (10) 

where D° represents the previous estimate in Eq. (jSj). It is apparent that the contribution 
of the polarization effect is negligible for temperatures T ^ a}. However, at temperatures 
comparable to aj, the contribution could be significant; notably, at those temperatures Fermi 
effects also need to play a role. 
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The isospin diffusion induced by mechanical forces has analogy in an electric current 
induced by the electric fields. Indeed, for large enough systems, the Coulomb interactions 
can contribute currents altering the concentration and, for completeness, we evaluate the 
conductivity aE for nuclear matter, relating the isospin flux to the electric field, 

Ti = aES, (11) 

where S is the local electric field. 



III. FLUXES FROM THE BOLTZMANN EQUATION SET 

A. Coupled Boltzmann Equations 

The two components of the binary system will be described in terms of the quasiparticle 
distribution functions /j(p, r, t). The local macroscopic quantities h{r,t) are expressed as 
momentum integrals of /, 

h{r, t) = I dp x(p)/(p, r, t) , (12) 

where g is the intrinsic degeneracy factor. Different standard expressions for macroscopic 
quantities in terms of /, such as pressure and heat flow, are listed in the Appendix 1X1 
The components are assumed to follow the set of coupled fermion Boltzmann equations, 

f + + F,-fi = J.. (13) 

ot rrii or op 

The terms on the l.h.s. account for the changes in /j due to the movement of quasiparticles 
and their acceleration under the influence of mean-field and external forces, included in Fj, 
while the r.h.s. accounts for the changes in /j due to collisions. In the following, we shall 
often denote the l.h.s. of a Boltzmann equation as Vi. With da/dQ and v* representing the 
differential cross section and relative velocity, respectively, the collision integral for particle 
1 is 



Jl — Jll + Jl2 



2(27rn)3 
, 9 



I d'p,^ dn' V* {J\ ha f[ fL - fl ha K Ka) 

- 1 d'p2 dn' V* (^^^ (A h f[ f, - h f2 f[ h) . (14) 
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Here, / = 1 — / is the Pauli principle factor. The factor of 1/2 in front of the first r.h.s. 
Jii term, compared to the J12 term, compensates for the double-counting of final states 
when integration is done over the full spherical angle in scattering of identical particles. The 
subscript a and the primes in combination with the particle subscripts 1 and 2 are used to 
keep track of incoming and outgoing particles for a collision. Other than in the context of 
particle components, such as here, the 1 and 2 subscripts will not be utilized in the paper. 
The collision integral J2 for particles 2 follows from ()14|1 upon interchange of the indices 
1 and 2. As it stands, the set of the Boltzmann equations (fT^. with (fT^ . preserves the 
number of each species. 

In the macroscopic quantities |T2|l . the distribution function / gets multiplied by the 
degeneracy factor g. When considering changes of macroscopic quantities (fT^ dictated by 
the Boltzmann equation (jlHj) . the changing distribution function / continues to be multiplied 
by (7. In the equation, the factor of / for the other particle in the collision integral J is 
accompanied by its own factor of g. As a consequence, in the variety of physical quantities 
we derive, the factor of / is always accompanied by the factor of (?, while, however, / is not. 
To simplify the notation, in the derivations that follow, we suppress the factors of g, only 
to restore those factors towards the end of the derivations. 

When the Boltzmann equation set is used to study the temporal changes of densities of 
the quantities conserved in collisions, i.e. number of species, energy and momentum, local 
conservation laws follow. Those conservation laws are discussed in Appendix |B] 



B. Strategy for Solving the Boltzmann Equation Set 

Irreversible transport takes place when the system is brought out of equilibrium such as 
in effect of external constraints. Aiming at the transport coefficients, we shall assume that 
the deviations from the equilibrium are small, of the order of some parameter e that sets the 
scale for temporal and spatial changes in the system. Then the distribution functions may 



be expanded in the power series in e 
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/ = + + f'^ + . . . , (15) 

where /'•'^^ represent the consecutive terms of expansion and /'•°'* is the strict local equilibrium 
solution. The terms of expansion in / may be nominally found by expanding the collision 
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integrals in e, following (fT3jl . expanding, simultaneously, the derivative terms in the equations 
and by demanding a consistency, 

Vf'+vf^ + ... = jf^ + j^ + jf^ + .... (16) 

Here, we recognize that the derivatives, themselves, bring in a power of e into the equations 
and, thus, the derivative series starts with a first order term in e. 

While we nominally included the zeroth-order term in the expansion of the collision 
integral Jj, the integral vanishes for the equilibrium functions 



= \ (17) 



A'j 



exp ^ — ^-ij, — y ~^ 

where /Xj, v and T are the local kinetic chemical potential, velocity and temperature which 
are functions of r and t, consistent with the Euler equations ()B7|) . Notably, the vanishing of 
the collision integrals is frequently exploited in deriving the form of the equilibrium functions, 
leading to the requirement that fj/fj is given by the exponential of a linear combination 
of the conserved quantities. In the context of specific transport coefficients, the boundary 
conditions for the Euler equations ()B7|) may be chosen to generate just those irreversible 
fluxes, and forces driving those fluxes, that are of interest. 

The equation set (fT^ can be solved by iteration, order by order in e, requiring 

= jf ) . (18) 

Thus, may be introduced into P,, producing and allowing to find f^^\ Next, 
inserting /j^^ into T>i yields T>f''' that allows to find /j^'' and so on. 

For finding the coefficients of linear transport, only one iteration above is necessary, since 
f^^\ as linear in gradients, yield dissipative fiuxes that are linear in those gradients. The 
local equilibrium functions on its own produce no dissipative fiuxes, as the species local 
velocities Vj and heat fiux Q vanish, while the kinetic pressure tensor P is diagonal, 

".Vf = /7^^/f(p,r,t)=0, (19a) 



(27r^)3 m 
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in the frame where the local velocity vanishes v(r, t) = 0, with E_ representing the local 
kinetic energy per particle. The above fluxes reduce the local continuity equations to the 
ideal-fluid Euler equations. 



C. Boltzmann Set in the Linear Approximation 

We now consider the terms linear in derivatives around a given point, i.e. the case of 
A; = 1 in ()18|) . for the Boltzmann equation set. On representing the distribution functions as 
fj = /j^"* + fj^\ we expand the collision integrals Ji, to get terms jf^'' linear in fj^\ Upon 
representing /j^"* as fj^^ = /j"'' (pj, we get for the k = 1 i = 1 version of (fTHj): 

^ + -^-^ + Fi-^ = -/n(0) - /i2(0) , (20) 
at nil or op 

where 

1 f '^^Pj'^ ^n' r.,* ( '^^v \ AO) AO) m' m' 



and where we have utilized the property of the equilibrium functions 

AO) AO) A0)> AO)' ^ AO) Ao) Aoy M)' /22^ 

J i J ja J i Jja J i J ja J i Jja' \ ) 

The result for i = 2, analogous to ()20|). is obtained through an interchange of the indices 1 
and 2. 

The l.h.s. of Eq. (l2Up contains the derivatives of equilibrium distribution functions with 
respect to t, r and p. These derivatives can be expressed in terms of the parameters describ- 
ing the functions ()17p . i.e. /i,, T and v. Through the use of the Euler equations (Appendix 
IB]) and equilibrium identities (Appendix [n)), moreover, the temporal derivatives may be 
eliminated to yield for the rescaled l.h.s. of (PUI) 



T / 5 \ p dT pp d p , 



\2mi 3 J rriiT dr mi dr pi 
Here, a symmetrized traceless tensor is defined as 5ry = | (xy + y x) — | (x • y) 1, and 
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where Sj is the entropy per particle for species i, = {5E_J3 — fii)/T. The result for species 
2 in the Boltzmann equation is obtained by interchanging the indices 1 and 2 in Eqs. 
and (1211). Note that dai = -di2. 

The representation ()23j) for the l.h.s. of the linearized Boltzmann equation ()2()j) exhibits 
the thermodynamic forces driving the dissipative transport in a medium. Thus, we have the 



tensor of velocity gradients contracted in with the tensor from particle momentum. 
The distortion of the momentum distribution associated with the velocity gradients gives 
rise to the tensorial dissipative momentum flux in a medium. As to the vectorial driving 
forces, they all couple to the momentum in 1^!^ and they all can contribute to the vector 
fluxes in the medium, i.e. the particle and heat fluxes, as permitted by the Curie law. 
The criterion that we, however, employed in separating the driving vectors forces in (j23p 
was that of symmetry under the particle interchange. When considering the diffusion in a 
binary system, with the two components flowing in opposite directions in a local frame, one 
expects the driving force to be of an opposite sign on the two species. On the other hand, 
in the case of the heat conduction, one expects the driving force to distort the distributions 
of the two species in a similar way in the same direction. 

Regarding the antisymmetric driving force in ()24|) . we may note that for conservative 
forces we have 

F. = -|-f/.. (25) 



We can combine then the first with the second term on the r.h.s. of ()24|1 by introducing the 
net chemical potentials /i* = /Xj + Ui and getting 



, PlP2 
ll2 - 



\- 








--) 


dT 


dr 




m2/ 






dr 



(26) 



P 

For a constant temperature T, the driving force behind diffusion is the gradient of difference 
between the chemical potentials per unit mass, /i*2 = p\/fni — ^\/m2., as expected from 
phenomenological considerations [igI. However, the temperature gradient can contribute 
to the diffusion as well, which is known as the thermal diffusion or Soret effect. We note 
that the vector driving forces in vanish when the temperature and the difference of net 
chemical potentials per mass are uniform throughout a system. 

Given the typical constraints on a system, it can be more convenient to obtain the driving 
forces in terms of the net pressure P*, temperature T and concentration 5, rather than /i^2 
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and T. Thus, on expressing the potential difference as /i5_2 — /^i2(-P*) ^) S^^ 



and 

J P1P2 
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(nfa VP* + n?2 VT + V5) , (28) 



p 

that we will utilize further on. The coefficient functions are 

and specific expressions for those functions in the nuclear-matter case are given in Appendix 
IdI Notably, however, the concentration 6 may not be a convenient variable in the phase 
transition region where the transformation between the chemical potential difference and 5 
is generally not invertible. 

With the l.h.s. of the linearized Boltzmann set (j2m) linear in the driving forces exhibited 
on the r.h.s. of (j^Hj) . and with the collision integrals linear in the deviation form-factors (j), 
the form factors need to be linear in the driving forces, 

4>i = -Ai ■ VT - fii : W - Ci ■ di2 , 

02 = -A2 ■ VT - I2 : W - C2 ■ di2 , (30) 

where A, B and C do not depend on the forces. On inserting into (pUjl. we get the 
following equations, when keeping alternatively a selected exclusive driving force finite: 

-^/r/r = /n(C) + /12(C), (31a) 
PiT 

^ A'^fP = /22(C) + /21(C), (31b) 



when keeping di2. 



P2T 



-P^/f)/T)=/n(/?) + /i2(/?), (32) 

nil 



and another one, with indices 1 and 2 interchanged, when keeping Vv, and, finally, 
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and another one, with 1 and 2 interchanged, when keeping VT (while di2 = 0). 

The hnearized coUision integrals lij cannot change the tensorial character of objects upon 
which they operate. Moreover, the only vector that can be locally utilized in the object 
construction is the momentum p. This implies, then, the following representation within 
the set dSni): 





= Q(p2) 


p 


(34a) 




pi ' 




A, 




( P' 5 \ p 

\2mi 3-7 miT2 ' 


(34b) 


1 


= hip') 




(34c) 



Here, the tensorial factors are enforced by construction. The factorization of the scalar 
factors is either suggested by the respective linearized Boltzmann equation or serves con- 
venience later on. The unknown functions a, b and c can be principally found by inserting 
into Eqs. The resulting equations are, however, generally quite complicated 

and analytic solutions are only known in some special cases. In practical calculations, we 
shall contend ourselves with a power expansion for the unknown functions. It has been 
shown that any termination of the expansion will produce lower bounds for the transport 
coefficients and that the lowest terms yield a predominant contribution to the coefficients 

HI- 



D. Formal Results for Transport Coefficients 

Before solving Eqs. ()31|) - (jnnj) . we shall obtain formal results for the transport coefficients, 
assuming that solutions to (j^ - ljH^ exist. We shall start with the diffusion. The velocity 
for species 1 is 

_ I f d'^P P _ 1 f d^P P , .(0) ir(O) 

Til J [znhy^ mi rii J [zirn)-^ mi 

-di2 I / C, . [/„(C) + In{C)] , (35) 
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where we have utihzed ()3U|) and ()34|) . The contribution of a tensorial driving force to 
the vector flow drops out under the integration over momentum, as required by the Curie 
principle. With a result for analogous to we get for the difference of average 

velocities (utilized for the sake of particular symmetry between the components) 

Y,-v,^-vt|{/^a,.mc) + /.(C)] 

"^'^ A2- [/22(C) + /21(C)] 

"^'^ Ci-[/n(C) + /12(C)] 




= -|({A,C}VT+{C,C}di2) , (36) 

where the brace product {■, ■} is an abbreviation for the integral combinations of vectors A 
and C, mult iply ing the driving forces. The brace product has been first introduced for a 
classical gas jll]. The fermio n g eneralization of the product and its properties are discussed 
in the Appendix |E| see also |l4| . 

The diffusion coefficient is best defined with regard to the most common conditions under 
which the diffusion might occur, i.e. at uniform pressure and temperature, but varying 
concentration. We have then, cf. Q, 

V,-V, = — rs = -!^D,VS, (37) 

(mi + 1712) ni n2 pi p2 

where mi2 is the reduced mass, l/mi2 = 1/mi + l/m2. Respectively, when P* and T vary, 
with di2 given by (j28|l . we write the r.h.s. of ()36|1 as 

Y, - V, = Ds (V6 + ^ VP* + kr Vt) , (38) 

P1P2 V / 

where, simplifying the notation, we dropped the subscripts 12 on coefficients 11. The diffusion 
coefficient in the above is given by 

3mi2 n \ p J 

and 

1 p {A,C} 
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We can note that the expressions above contain 11'' in denominators. Normally, positive 
nature of the derivative ()29cp is ensured by the demand of the system stability. However, 
across the region of a phase transition the concentration generally changes while the chemical 
potentials generally do not, so that 11'' = 0. While the coefficient above is the one we 
are after as the standard one in describing diffusion, in the phase transition region it can 
be beneficial to resort to the description of diffusion as responding to the gradient of the 
potential difference in (j2HI)- Notably, as explained in the Appendix 10 the brace product 
{C, C} in ()39p is positive definite. This ensures the positive nature of away from the 
phase transition and, in general, ensures that, at a constant temperature, the irreversible 
asymmetry fiux fiows in the direction from a higher potential difference /i*2 to lower. 

As to the Soret effect, i.e. diffusion driven by the temperature gradient, described in 
fl38|) - (jim) . it has its counterpart in the heat flow driven by a concentration gradient, termed 
Dufour effect. Transport coefficients for counterpart effects are related through Onsager 
relations l3] that are also borne out by our results. The diffusion driven by pressure is 
rarely of interest, because of the usually short times for reaching the mechanical equilibrium 
in a system, compared to the equilibrium with respect to temperature or concentration. 
However, an irreversible particle flux may be further driven by external forces, such as due 
to an electric fleld £. With the flux induced by the fleld given by Vs = cte^, where cte 
is conductivity, with the flrst equality in ()37|) . and with (jHUj) and (j21I), we flnd for the 
conductivity 

..^-^fElPlY f^.SL) {C.C}= (^-A) " B„ (41) 

3mi2 \ P J \"^2 ^1 J \"^2 rrii J ff 

where is charge of species i. We see that conductivity is closely tied to diffusivity. 

While our primary aim is to obtain coefficients characterizing the dissipative particle 
transport, due to the generahty of our results we can also obtain the coefficients for the 
transport of energy and momentum. Thus, starting with the expression ()Ale|) in a local 
frame and proceeding as in the case of ()35p and ()36p. we get, with ()33|). 

Qi + Q2 = -| ({A, A} VT + {C, A} di2) + ^ (^1 ni Vi - V2) 

= -| ({A, A} VT + {C, A} di2) (42) 

^(E, ^2^ P1P2 



3 V rrii m2 
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where in the second step we make use of the condition on local velocities pi Vi + P2^2 = 0. 
The standard procedure in coping with the heat flux is to break it into a contribution 
that can be associated with the net movement of particles and into a remnant, driven by 
the temperature gradient, representing the heat conduction. With this, the driving force 
di2 needs to be eliminated from the heat flux in favor of the species velocities. Using (jHEI), 
we find 



Q1 + Q2 = -VT^f{A,A} 



3 V {C,C} 
+ (Vi - V2 



5 [E, E,\ pip2^{C,A} 



mi 1712/ p {C, C} 



(43) 



..^f{A,A}-i?^), (44) 



The coefficient 

" " 3 \^L-'-J {C,C} 
relating the heat flow to the temperature gradient, is the heat conduction coefficient. From 
(j44p and considerations in Appendix 10 it follows that k, given by Eq. (j44j) is positive definite. 

The final important coefficient that we will obtain, for completeness, is the viscosity. The 
modification of the momentum flux tensor ()Ald|) . on account of the distortion of momentum 
distributions described by (IHUj). is 

=(1) ^ r I [ _fP_^xf 

J {2T^hf mi J {2Tihf m2 

d^^P ¥^ H AO) m f _d^Wf (= =\ (0) y(o) 



{2'Khf mi V y •'i •'1 J [2'Khf ma 

f f ■ ^ ^ f . ^ f(0) f(0) 



= -|vv{i?,i?}. (45) 

The coefficient of proportionality between the shear correction to the pressure tensor and 
the tensor of velocity derivatives is, up to a factor of 2, the shear viscosity coefficient 

77 = ^{1,1}. (46) 

As with other results for coefficients, from Appendix lEl it follows that the result for rj above 
is positive definite, as physically required 

On account of symmetry considerations within the linear theory, the changes in temper- 
ature or concentration do not affect the pressure tensor. However, the situation changes if 
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one goes beyond the linear approximation. For a general discussion of different higher-order 



effects see Ref. 



ll|. As a next step, we need to find the form factors in ()30|): that requires 



finding the functions a, b and c in by solving Eqs. (|!?T|) - (p?!^ . 



IV. TRANSPORT COEFFICIENTS IN TERMS OF CROSS SECTIONS 

A. Constraints on Deviations from Equilibrium 

Since the zeroth-order, in derivative expansion, local-equilibrium distributions are con- 
structed to produce the local particle densities, net velocity and net energy, corrections 
to the distributions cannot alter those macroscopic quantities. Thus, we have locally the 
constraints 

= I (^^"'-°- 

With driving forces being independent of each other and with form factors in (jHUj) being in- 
dependent of the forces, each of the form factors sets must separately meet the constraints. 
By inspection, however, one can see that the density and energy constraints are met au- 
tomatically with the forms (|34p of form factors. Moreover, the tensorial distortion (|34cp 
satisfies all the constraints. At a general level, the ability to meet the constraints while solv- 
ing Eqs. (jHT|) - (jnn|) relies on the fact that the linearized collision integrals lij (in Eqs. ()2Up and 
(j2^) nullify quantities conserved in collisions, so a combination of the conserved quantities 
may be employed in constructing the form factors (pi, ensuring that the constraints are met. 
When the transport coefficients get expressed in terms of the brace products, though, the 
ensuring that the constraints are met becomes actually irrelevant for results on the transport 
coefficients, because the linearized integrals and correspondingly brace products nullify the 
conserved quantities. 

Given cross sections and equilibrium particle distributions, the set of equations (jHT|) - (jn^ 
may be principally solved. However, such a solution is generally complicated and would 
likely not produce clear links between the outcome and input to the calculations. On the 
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other hand, the experience has been that when expanding the form-factor functions, a, b and 
c in (jnH), in power series in p^, the lowest-order results represent excellent approximations 
to the complete results and are quite transparent, e.g. Thus, we adopt here the latter 
strategy and test the accuracy of our results in a few selected cases. 



B. DifFusivity 

If we insert ()34a|) with Cj(p^) = q into the local velocity constraint ()47ap . we get the 
requirement 

^ ■ /f ^ ^ + - / 77^ ir fT = . (48) 



After partial integrations, we find that this is equivalent to the requirement Ci = — C2 = c. 

When Ci is constant within each species, then Cj is up to a factor equal to momentum 
and, thus, gets nullified by the linearized collision integral within each species /m(C) = 0. 
To obtain a value for c, we multiply the first of Eqs. ^T\i by Ci and the second by C2, 
add the equations side by side and integrate over momenta. With this, we get an equation 
where both sides are explicitly positive definite and, in particular, the l.h.s. is similar to the 
l.h.s. of Eq. ()48|) . but with an opposite sign between the component terms. That side of 
the equation can be integrated out employing the explicit form of f^^^ from Eq. (|T7jl . The 
other side of the resulting equation represents {C, C} where only the interspecies integrals 
survive. On solving the equation for c, we find 

6piP2 



PX12 

where 



(49) 



2 / d^pi d^p2 ^(dai2\ I /n2 /(o) .(0) ?(o)/ ?(o)/ 



The integral stems from a transformed brace product {C,C} and we resurrect here the 
degeneracy factors g. For the brace product itself, we find 

c^ ( p V 18 



{C,C} = - \^\ xi2 = —. (51) 
2 \P1P2J X12 



On inserting this into the diffusivity (jSHl), "we obtain 

6T n-^ /piP2^' 



Ds = ^^ ^ . (52) 

mi2 nxi2 \ P ' 
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In the above, we see that the diffusion coefficient depends both on the equation of state, 
through the factor 11'', and on the cross section for colhsions between the species, through 
Xi2- The collisions between the species are weighted with the momentum transfer squared. 
Only those collisions between species that are characterized by large momentum transfers 
suppress the diffusivity and help localize the species. The marginalization of collisions with 
low momentum transfers is a common feature of all transport coefficients. 

At high temperatures the Fermi gas reduces to the Boltzmann gas. In absence of mean- 
field effects, we find 11'^ ~ — for small asymmetries. The integral Xi2 is then of the order 



n'^cri2p^/m ~ r? a \fmT^ . Together, these yield -D5 ~ ^T^y m' '^'^^ precise hig i-T result 
for isotropic cross-sections in the interaction of species with equal mass m is (ill . Il5l | 



D, = -^^\—. (53) 

The square-root dependence on temperature will be evident in our numerical results at high 
T . With an inclusion of the mean field, with the net energy quadratic in asymmetry, the 

derivative 11'^ gets modified into II'^ ~ -^^ —. Thus, the mean field enhances the diffusion. 

At low temperatures, the derivative 11'' is simply proportional to the symmetry energy, 
11'' ~ As to the collisional denominator of the diffusion coefficient, at low temperatures 
the collisions take place only in the immediate vicinity of the Fermi surface. We can write 
the product of equilibrium functions in the collision integral as 



fTfffT'fT' = K,K,K',K',, where K, = , (54) 



1 



2 cosh ^^^^^ — 



and, at low T, Ki ~ 211 vnT bij?' — "p^p^. The integration in (j5U|) yields X12 ~ oyi"^ T"^ lv\ 

2 



In consequence, we find that the diffusion coefficient diverges as 1/T at low temperatures. 
For the spin diffusion coefficient, one finds within the low temperature Landau Fermi-liquid 



theory 



1^ 



= y (1 + K) , (55) 



where fi? is Fermi velocity, Fq is a spin-antisymmetric Landau coefficient and To is a char- 
acteristic relaxation time that scales as To ~ T~^. The isospin diffusivity for symmetric 
matter should differ from the spin diffusivity in the replacement of the spin-antisymmetric 
Landau parameter with the isospin asymmetric parameter, neither of which has a significant 
temperature dependence. Thus, here consistently we find a divergence of the diffusivity 
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at low temperatures. Moreover, the factor (1 + F^) is nothing else but a rescaled symmetry 
energy, with Fq being the ratio of the interaction to the kinetic contribution to the energy 
jlfll . Thus, here consistently we find a proportionality of the diffusivity to the symmetry 
energy at low temperatures. 

To summarize the above results on diffusivity, we find that the diffusivity is inversely 
proportional to the cross section between species for high momentum transfers. Moreover, 
whether at low or high temperatures, the diffusivity is sensitive to the symmetry energy 
in the mean-fields. The mean-field sensitivity is associated with the factor 11'' = + 
^ (ml ~ ^) ~ + ■where the last equality pertains to the system of neutrons 
and protons and represents the interaction contribution to the symmetry energy at the 
relevant density. 

While we obtained the diffusivity here assuming constant Cj in ()34ap . we will show that 
the next-order term in the expansion of Cj increases the diffusion coefficient Ds only by 2% 
or less in our case of interest. 



C. Heat Conductivity 

Evaluation of the heat conduction and shear viscosity coefficients requires similar method- 
ology to that utilized for the diffusivity. While these coefficients have been obtained in the 
past for a one component Fermi system Jt^, it can be still important to find them for the 
two component system. 

If we assume aj(p^) = in fl34b|l . then, interestingly, we find that the momentum con- 
straint ()47a|) is automatically satisfied. To obtain the values for Oj, we multiply Eq. ()33p on 
both sides by Ai and integrate over momenta and we multiply the equation analogous to 
dnSl) by A2 and also integrate over momenta. As a consequence, we get a set of equations 
for of the form 

Lj = Ajiai + Aj2a2 , J = 1,2, (56) 
where Aji are coefficients independent of a, 

A^ = \ ([A, A]ii + [Ai, AJ12) , A12 = A21 = ^— [Ai, A2]i2 , (57) 

ttj ai a2 

cf. Appendix IE} and 

1 / _o 25 ._.9\ 
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where Ej and {E_jY are, respectively, the average local square kinetic energy of species j 
and square average local kinetic energy of the species. 
The solution to the set (jHHjl is 

ai = {A22 Li - A12 L2) /A^ , 

a2 = {AnL2-Ai2Li)/A^, (59) 

where the determinant is 

Aa = AuA22- AI2. (60) 
The brace product {A, A} for use in calculating the heat conduction coefficient k, in ()44j) is 

{A, A} = aiLi + a2i^2. (61) 

The product {C, A} in (jUj) can be calculated given the values of a and c, and {C, C} was 
already obtained before. 

D. Shear Viscosity 

Evaluation of the shear viscosity coefficient t] follows similar steps to those involved in 
the evaluation of k. Thus, we assume hi{p^) = hi in (|34c|) . To find the coefficient values, we 
convolute both sides of Eq. ()H2|) with p p and integrate over the momenta and we do the 
same with the other constraint equation for B. The l.h.s. integrations produce 

d^P = = AO) r(o) 2 f d^P 4 AO) m 20 

PP : PP fi fr = o ^ / TTT^P fi fi =^PiEi. (62) 



TTiiT J (27rn)3*^*^ • 3miT J 3 

With the above, we get the set of equations for hi. 
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— PjE^=B,^h^ + B,2h2, J = l,2, (63) 



where the coefficients E are given by. 



^ii = \PP, PP\ii + [(P P)i, (P P)i]l2 , B12 = B21 = [(p p)l, (p p)2]l2 • (64) 

Solving the set for h, we find 

20 

bi = ■:ri^{piE^B22- P2E2B12) , 
20 

h = ^ (P2 ^2 Bu - pi E^ S12) , (65) 
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where the determinant is 

AB = i3ni322-SL- (66) 
The brace product for calculating the shear viscosity coefficient r] = j^{B,B} becomes 

= = 20 

{B,B} = - (6i pi E, + b2 P2 E,) . (67) 

V. QUANTITATIVE RESULTS 
A. Transport Coefficients 

We next calculate the transport coefficients as a function of density and temperature, 
using experimentally measured nucleon-nucleon cross sections. The cross sections may be 
altered in matter, compared to free space, but the modifications are presumably more im- 
portant at low than at the high momentum transfers important for the transport coefficients. 
With regard to the diffusivity, we first ignore any mean-field contribution to the chemical po- 
tential difference between species. This yields a reference diffusivity to which the diffusivity 
affected by mean fields may be compared. 

The diffusivity for the experimental cross sections and no interaction contributions to 
the symmetry energy is shown at 5 = and different densities n in Fig. Q as a function 
of temperature T. At low temperatures, the diffusivity diverges due to a suppression of 
collisions by the Pauli principle. At high temperatures, compared to the Fermi energy, the 
role of the Pauli principle is diminished and the diffusivity acquires a characteristic ^/T 
dependence. At moderate temperatures and densities in the vicinity and above normal, the 
diffusion coefficient turns out to be in the vicinity of our original estimate of Dj ~ 0.2 fmc. 

It should be mentioned that, for symmetric matter, the factors for temperature and 
pressure gradients in the thermodynamic force di2 (I28p vanish, 11'^ = and 11"^ = 0, and the 
brace product in (jiUI) vanishes, {A, C} = 0, yielding A;^ = in (jHHj) . As physically required, 
the temperature and pressure gradients produce no relative motion of neutrons and protons 
for the symmetric matter. 

The diffusivity at normal density at different asymmetries is next shown in Fig. |21 as a 
function of temperature. Because of charge symmetry, the diffusivity does not depend on 
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the sign of S. At low temperatures the diffusivity is generally expected to behave as 

Di 4lF ' (68) 
while at high temperatures in the manner prescribed by ()53p. With the respective behaviors 
serving as a guidance, we provide a parametrization of our numerical results for Ds as a 
function of n, T and 6, 

1.54 



Dr = (1-0.195^1 



11.34 / n 



T2-38 \no 



1.746 f n V'^^ , r. 0.^0.913 f ^0 



+ — + 0.00585 T^-^^'^ — 

T \noJ \n 



(69) 



Here, the temperature T is in MeV and the diffusivity Dj is in fmc. The parametrization 
describes the numerical results to an accuracy better than 4% within the region of thermo- 
dynamic parameters of 1.0 < n/no < 4.0, 10 MeV < T < 100 MeV and \5\ < 0.4. This is, 
generally, the parameter region of interest in intermediate-energy reactions. 

The heat conductivity is shown for symmetric matter at different densities in Fig. El as 
a function of temperature. The results are similar to those in Ref. |3], though there the two 
component nature of nuclear matter was ignored and the isospin-averaged nucleon-nucleon 
cross-sections have been used. A closer examination of results in Subsections II V (^1 and II V Dl 
indicates that the use of the isospin-averaged cross-sections is, actually, justified for symmet- 
ric matter, when calculating the heat-conduction and shear-viscosity coefficients. Otherwise, 
however. Fig. El has been based on a more complete set of cross sections than results in 
As in the case of diffusivity, the heat conductivity diverges at low temperatures and tends 
to a classical behavior at high temperatures, exhibiting there no density dependence and 
being proportional to velocity, k oc VT. As in the case of diffusivity, we next provide a 
parametrization of our numerical results for the heat conductivity k as a function of n, T 
and 6, 



:i + 0.105^ 



0.235 / n - 



(70) 



0.0816 / \ 0.0171 

-0.0582 ( — I + 0.0238T°-^^2^ ( — 

Here, T is again in MeV and k, is in c/fm^. The parameterization agrees with the numer- 
ical results to an accuracy better than 4% within the range of thermodynamic parameters 
indicated in the case of Dj. 
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The shear viscosity coefficient t] is shown for symmetric matter at at different densities, as 
a function of temperature, in FigHJ Again, the results are similar to those in Ref. f^. At high 
temperatures, the dependence on density weakens and the viscosity becomes proportional 
to velocity. The numerical results for t] are well described, to an accuracy better than 4% 
within the before-mentioned range, by 



7] = (1 + 0.105^ 



856 / n V-^^ 240.9 / n ^ ^'^^ 



2^1542-0.76 



(71) 



Here, rj is in MeV/fm^ c and T is in MeV. 

We note in ()69| ) - ()7H) . that the diffusion coefficient weakly drops with increasing magnitude 
of asymmetry \6\, while the viscosity and heat conduction coefficients weakly increase. Given 
the weaknesses of the dependencies, the behaviors exhibited in parametrizations represent, in 
practice, averages over the considered independent-parameter regions. Overall, the drop and 
rise in the respective coefficients with |5| is characteristic for a situation where the local flux of 
a component grows faster than the concentration of that component. That type of growth, 
with the magnitude of asymmetry, typifies a mixture of degenerate fermion gases. The 
general trends can be deduced following the mean-free-path arguments from Sec. |nj When 
the average velocity rises with asymmetry, so do the heat conduction and shear viscosity 
coefficients. Additional rise for those coefficients, in the case at hand, can result from the 
Pauli principle effects and from the difference between cross sections for like and unlike 
particles. Regarding the diffusion coefficient, though, one needs to consider an irreversible 
part of relative particle flux, under the condition of the concentration varying with position. 
If, starting with a given configuration of concentration gradients, one introduces uniform 
changes of concentration on top, not just the overall relative flux undergoes change but also 
the reversible flux of concentration gets altered. The rise in the relative flux associated 
with the velocity of a dominant component rising with concentration is normally more than 
compensated by the rise in reversible flux, leading to a reduction in the irreversible flux and 
producing a reduction in diffusivity with particle asymmetry. A mean-field example where 
the reversible flux eats into the net flux reducing the diffusivity with increasing asymmetry 
is the estimate in Eq. (jlOj) . obtained there without invoking the particle Fermi statistics. 

As is found in Sees. IIII (^1 and IIVB| the dependence of mean fields on species enters 
the diffusivity through the factor 11'' resulting from the variable change in thermodynamic 
driving force, from the difference of chemical potentials per mass to asymmetry. The simplest 
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case where one can consider the impact of the mean fields is that of the symmetric nuclear 
matter, at 5 = 0. In this case, the factor may be represented as 



where = ^ = ^ (cf. Appendix E)). At high temperatures, we have approximately 
so that I ~ 2T. The naive expectation is that has a linear dependence on the net 
density, = of' ? where af = 14 MeV and u = 1. The mean-field amplification factor 

R = n^(a^)/n^(a^ = 0) for the diffusion coefficient, assuming the linear and also quadratic 
density-dependence of (z/ = 1 and 2) is shown in Fig. El The quadratic dependence gives 
higher amplification factors at n > hq, than the linear dependence, while the opposite is 
true at n < uq. At low temperatures and moderate to high densities the amplification is 
very strong suggesting that the diffusion could be used to probe the symmetry energy, aside 
from the in-medium neutron-proton cross sections. 



The calculations of transport coefficients above have been done assuming that the func- 
tions Oj, bi and q in Eqs. (|34j) can be approximated by constants. In the more general case, 
the functions can be expanded in the series in p"^, e.g. 



The coefficients of the expansion can be found by considering moments of the form-factor 
equations (jHT|) - ()33|) . With the more general form of the form-factor functions, the transport 
coefficients generally increase, but their rise is generally very limited. 

To illustrate the magnitude of higher-order effects, we provide in Table H] results for the 
diffusivity obtained in the standard first-order and in the higher-order calculations at sample 
densities and temperatures. In the indicated cases, the second-order calculations never 
increase the diffusion coefficient by more than 3% above the first-order calculations. The 
efficiency of our Monte-Carlo procedure employed to evaluate the integrals for coefficients 
worsens as the order of our calculations increases and, correspondingly, we provide only a 
single third-order result for illustration. 




(72) 



B. Testing the Form-Factor Expansion 




(73) 
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C. Isospin Equilibration 

To test sensibility of the results on diffusivity and to gain an elementary insight into the 
process of isospin equilibration in a reaction, we carry out a schematic consideration of the 
equilibration. For definiteness, and to ensure a level of applicability for our consideration, 
we take the case of a ^^Ra + ^^Zr reaction at Eiab/A = 100 MeV. Densities in the central 
region of the reaction are not far from normal. Following the degenerate Fermi-gas limit, 
the temperature in the central region can be estimated with T ~ a/ -E'zab/ (2a) ~ 20 MeV 
using a ^ A/ {8 MeV); deviations from the degenerate limit yield a bit higher value. Under 
those conditions, basing on Figs. Q and El we estimate the streaming contribution to the 
diffusivity in the central region at 0.21 fmc and the mean-field contribution at 0.20 fmc, for 
a net Dj ^ 0.41 fmc. 

Considering the direction perpendicular to the plane of contact between the nuclei, with 
nuclei extending a distance L ~ (A/noY^^ ~ 8 fm both ways from the interface, we may use 
the one-dimensional diffusion equation to estimate the isospin equilibration 

where x is the direction perpendicular to the interface, cf. Eq. (|7j). With isospin flux van- 
ishing at the boundaries of the region [— L, L], the solution to (fflj) is 

oo oo 

6{x, t) = 6oo + ttn sin kn X exp (— -D/ t) + 6„ cos g„ x exp {—Dj t) (75) 

n=l n=l 

where knL = — |) tt and qnL = rni. The coefficients a„ and 6„ are determined by the 
initial conditions and, in the case in question, 6„ = 0. 

The different terms in the expansion ()75|) correspond to the different levels of detail in the 
distribution of concentration, as characterized by the different wavevectors. We see that the 
greater the detail the faster the information is erased, with the erasure rates proportional to 
wavevectors squared. The overall distribution tends towards 5oo as t — > oo. The late-stage 
approach to equilibrium is governed by the rate for the term with the lowest wavevector, 
i.e. ai. Defining the isospin equilibration time tn as one for which the original isospin 
asymmetry between the nuclei is reduced by half, we get from ()75|) an estimate for the 
reaction 

ln2 41n2L2 _ , 
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for the case above. When we carry out the full respective Boltzmann-equation simulations 
of the lOOMeV/nucleon ^^Ru + ^^Zr reactions, at the impact parameter of 6 = 5fm > L/2 
(to ensure a neck comparable to that in the consideration), we find that, indeed, the nuclei 
need to be in contact for about 40fm/c for the isospin asymmetry to drop to the half of 
original value. 

VI. SUMMARY 

Diffusion and other irreversible transport phenomena have been discussed for a binary 
Fermi system close to equilibrium. For weak nonuniformities, the irreversible fluxes are linear 
in the uniformities, with the characteristic transport proportionality-coefficients dependent 
only on the equilibrium system. It is hoped that, in an analogy to how the nuclear equation 
of state and symmetry energy are employed, the coefficient of diffusion can be employed to 
characterize reacting nuclear systems with respect to isospin transport. 

Following a qualitative discussion of the irreversible transport in the paper, the set of 
coupled Boltzmann-Uhlenbeck-Uehling equations was considered for a binary system, as- 
suming slow macroscopic temporal and spatial changes. The slow changes allow to solve the 
equation set by iteration, with the lowest-order solution being the local equilibrium distri- 
butions. In the next order, corrections to those distributions were obtained, linear in the 
thermodynamic driving forces associated with the system nonuniformities. These corrections 
produce irreversible fluxes linear in the forces. The transport coefficients have been formally 
expressed in terms of brace products of the responses of distribution functions to the driv- 
ing forces. The considered coefficients include diffusivity, conductivity, heat conduction and 
shear viscosity. 

The set of the linearized Boltzmann equations was, further, explicitly solved under the as- 
sumption of simplified distribution-function responses to the thermodynamic driving forces. 
The solutions to the equations led to explicit expressions for the transport coefficients, with 
the diffusivity given in terms of the collision integral for collisions between the two species 
weighted by the momentum transfer squared. Besides associated sensitivity to the cross sec- 
tion for collisions between the species, the diffusivity is also sensitive to the dependence of 
mean fields on the species. The collisions between the species are those that inhibit the rela- 
tive motion of the species; the difference between mean fields affects the relative acceleration 
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and, in combination with the coUisions, the stationary diffusive flux that is estabhshed. 

We calculated the isospin diffusivity for nuclear matter, using experimental nucleon- 
nucleon cross sections for species-independent mean-flelds. At low temperatures and high 
densities, the diffusivity diverges due a suppression of collisions by the Pauli principle. At 
high temperatures, the diffusivity is roughly proportional to the average velocity and is 
inversely proportional to the density The diffusivity weakly decreases with an increase in 
the absolute magnitude of asymmetry. We provided an analytic flt to our numerical results. 
For completeness, we also calculated the heat conduction and shear viscosity coefficients 
and provided flts to those. Moreover, we calculated the diffuseness mean-field enhancement 
factor for symmetric matter, assuming a couple of dependencies of the symmetry energy 
on density. At low temperatures, the enhancement factor is simply proportional to the 
net symmetry energy divided by the kinetic symmetry energy. Considering the expansion 
of the form-factors in distribution-function responses, we demonstrated that corrections to 
the Boltzmann-equation transport coefficients, beyond the approximations we employed, 
are small. Finally, we produced an elementary estimate for isospin equilibration in a low 
impact-parameter collision. 
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APPENDIX A: MACROSCOPIC QUANTITIES 

We shall consider different types of macroscopic quantities, either net or for separate 
components, either in the general frame of observation or in a local frame. For a single 
component i in the observation frame, the density rii, mean velocity Vj, mean kinetic energy 
Cj, momentum fiux tensor Pj and kinetic energy fiux qj, are given in terms of the distribution 
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fi, respectively, as 



ni{r,t) = j^^^ j d^PfiiP,r,t), (Ala) 
^|#,2.y.(p,,,<), (Alb) 



^2 



_^y,3pA.;,(p.,.,). (Ale) 



9 I PP 



(2^/<''f'^|-^-'P------ 

The net quantities result from combining the component contributions. Thus, the net 
density is n = rii + n2, the net mass density is p = pi + p2 = fni ni + 1712 n2 while the 
net velocity v is obtained from p v = pi + p2 The kinetic energy e averaged over all 
particles is given by n e = rii Ci + n2 63, the net momentum flux is p = Pi + P2 and the net 
kinetic energy flux is q = qi + q2. 

Local quantities are those calculated with momenta transformed to the local mass frame, 
i.e. following the substitution p ^ p — rrij v. To distinguish local quantities from those in 
the observation frame, when the frame matters, the local quantities will be capitalized. The 
local momentum flux tensor P is the kinetic pressure tensor and the local kinetic energy 
flux Q is the heat flux. 



APPENDIX B: CONTINUITY EQUATIONS 

The collisions in the Boltzmann equation set ()13p conserve the quasiparticle momentum 
and energy and the species identity. This leads to local conservation laws for the corre- 
sponding macroscopic quantities. 

Let Xi(p) represent one of the quasiparticle quantities conserved in collisions, Xj(p) = ^ij^ 
p or p^ /2mj. For those quantities, the integration with collision integrals produces 

J2 [ d'PXjJj = 0. (Bl) 
As a consequence, from the Boltzmann equation set, we obtain 



y 



x,ft + ^'t + Frf 1=0. (B2) 



^ J (27rn)3 \dt nij dr ^ dp 
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After a partial integration, we get from the above 



where the averages are defined with 

Substituting for Xj the conserved quantities {Xj{p) = ^iji P or P^/2^j), we get the 
respective continuity equations: 

|^ + |:-KvJ = (B5a) 

d d 

^(py) + — •p-niFi-n2F2 = 0, (B5b) 
ot or 

d d 

^ e) + — • q - ni Vi • Fi - 712 V2 • F2 = . (B5c) 
ot or 

Here, we made yet no use of the local frame. 

The local frame is useful when wants to make use of the assumption of local equilibrium 
that imposes restrictions on local quantities. On representing the average velocities as 
Vj = Vj + V in the equations above, we obtain the following set, 

^ + ^ • (n. v) + - • (n, VJ = , (B6a) 
| + |-(PX)^0, ^ (B6b) 

^ M + ■ (Pvv) + |- . P - niFi - n2F2 = 0, (B6c) 



g-^{nE) + — -{nEY) + P:—v + — -Cl-n,V^F^-n2V^F2^0. (B6d) 

The equation for mass density in the set above follows from combining the equations for 

particle densities. 

The above equations significantly simplify when the assumption of a strict local equi- 
librium is imposed. Under that assumption, the local species velocities and the heat flow 
vanish, Vj = and Q = 0, and the kinetic pressure tensor becomes diagonal, P — ^nE^l. 
The equations reduce then to the Euler set 

OT) ■ f) 

^ + ^-Kv) = 0, (B7a) 

d , , d , , 2d{nE) ^ ^ ^ x 

Qi (P^) + • ^) + 3 - rzi Fi - ^2 F2 = , (B7b) 

^(n^)+v-^(n^) + ^n^^-v = 0. (B7c) 
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APPENDIX C: SPACE-TIME DERIVATIVES FOR AN IDEAL FLUID 

In an ideal fluid, all local quantities can be expressed in terms of the local temperature 
T and the local kinetic chemical potential /ij. If we consider changes of the densities Ui 
or of the local kinetic energies with respect to a parameter x representing some spatial 
coordinate or time, or their combination, we find 

(-1 Si-' o "Tri^^O ' 

OX ox 2 ox 

where aj = /ij/T, /? = logT and = {drii/ dnijx- With the trace derivative defined as 

d _ d d 
dt dt dr ' 

a particular version of the above relations is 

drii dai 3 djd 

m ^ ~dt^2'''~di' 

dt 2 " dt 2 dt ^ ' 

A combination of the above trace-derivative relations with the Euler equations from 
Appendix IHl yields the following simple results, 

^ = , (C3a) 
2 d 

-dl = -3d-r-^' ^^'^^ 
the consistency of which with ()(y2|l and ()B7|1 is easy to verify. The results ()(y3j) express 
basic features of the isentropic ideal-fluid evolution of a mixture. The entropy per particle 
in species i depends only on a^, while the ratio of the densities of species rii/n^ depends 
both on ai and a2- The conservation of ai for both species is equivalent to the conservation 
of entropy per particle and of relative concentration. Finally, the density for species i is 
proportional to T^/^ multiplying a function of a^, which is equivalent to the second of the 
results above, given the continuity equation for species and the conservation of a,. 

APPENDIX D: VARIABLE TRANSFORMATION 

The driving forces for diffusion are naturally expressed in terms of the gradients of tem- 
perature and of chemical potential difference per unit mass /x^g- However, given the typical 
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constraints on systems, it can be convenient to express the chemical potential in terms of 
other quantities, that are easier to assess or control, such as the differential concentration 6, 
temperature T and net pressure P*. A transformation of the variables for the driving forces 
has been employed, at a formal level, in Sec. IIII CA Here, we show, though, how the trans- 
formation can be done in practice for the interaction energy per particle specified in terms 
of the particle density n and concentration 6, = E'"{n,5). With the nuclear application 
in mind, we limit ourselves to the case of mi = m2 = m. 

The transformation can exploit straightforward relations between different differentials. 
One of those to exploit is the Gibbs-Duhem relation 



dP^ = Til dfi\ + n2 + ns dT = n dfi'^ + 



mn5 



dfi\2 + ns dT . 



(Dl) 



Here, s is the entropy per particle and /i* = (yU* +/i2)/2 is the median chemical potential. Two 
other relations stem from the differentiations of equilibrium particle distributions, already 
utilized in Appendix El 



drii = djjLi + 



dT = ^idfii + 



dT. 



(D2) 



With fi^ = d{nE'")/dni, on adding and subtracting the two {i = 1, 2) relations side by side, 
we find 



dn 



djj: 



dn 



X 



t 

12 



9/i 



dn 



dn 



dn — 



~d5 



12 



85 



d6 



d6 



+ 



+ y (ei - 6) 



dni 



dT 



dn2 

w 



M2 



dT (D3) 



and 



6 dn + nd6 = — ^2) 



df/ - 



dn 



dn — 



dfx 



t 

12 



dn 



dn 



5 
86 



d6 



~d5 



+ 



d5 



dni 



AH 



dn2 
'df 



M2 



dT . (D4) 



Those two equations have the structure 



Gkn dn — Gkn dfi^ + Gkd 0^/^12 + Gks d6 + G^t dT , 



(D5) 



where k = 1,2 and where the coefficients G can be worked out from ()D3|) and ()D4|) . On 
multiplying the sides of the first {k = 1) equation by G2n and the sides of the second {k = 2) 
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equation by Gi„ and on subtracting the equations side by side, we can eliminate the dn 
differential obtaining 

= {G2n Gin ~ + {G2n G\d — Gin G2d) C^A*12 

+ (G^2n ^15 — Gin G25) dS + {G2n Gyr — Gin G2T) dT 

= R^dn^ + Rddn{^ + RsdS + RrdT . (D6) 
On eliminating next the diit differential using the Gibbs-Duhem relation, we find 



= ^ = . (D7c) 

APPENDIX E: BRACE ALGEBRA 

The brace products are employed in finding the transport coefficients within linear ap- 
proximation to the Boltzmann equation. The brace product of two scalar quantities A and 
B associated with the colliding particles is defined as 

= l>l,B],i + lAB]i2 + [ABL , (El) 

where, in the last step, we have broken the brace product into square-bracket products 
representing contributions from collisions within species 1, from collisions between species 1 
and 2 and from collisions within species 2, respectively. 

We will first show that the square-bracket product is symmetric. Thus, we have explicitly 

^^'^^^^ "27 (2^ (2^''^^'' ['dnj 

xAia {Bia + Bib — -Bj'a ~ -^ife) 

d^Pa d^pb f daii\ (0) (0) ?(o)/ ;(o)/ 

I ~lrr Jia Jib Jia Jib 



8 J (27r)3 (27r)3 V 
X ( Aa + A,b - - A'rb) (Bia + B^b - B^ - B'J , (E2) 
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where, to get the last result, we have first utilized an interchange of the particles in the 
initial state of a collision and then an interchange of the initial and final states within a 
collision. It is apparent that the r.h.s. of ()E2|1 is symmetric under the interchange of A and 
B. Moreover, we can see that a square bracket for B = A, [A, Aju, is nonnegative and that 
it vanishes only when A is conserved in collisions. 

We next consider the contribution from collisions between different species, 

\A m - ! d!^V2 rn,.* ( ^^12 A ,(0) .(0) 7(0)/ m, 

x{A^ + A2){B^ + B2-B[-B'^) 
_ 1 f d^pi d'^p2 ( dai2\ AO) AO) ?(o)/ ?(o)/ 

"27 j2^'Wr \dn) 

X {Ai + A2- A[ - A'^) {Bi + B2-B[- B'^) . (E3) 

Here, we again utilized an interchange between the initial and final states and we again 
observe a symmetry between A and B on the r.h.s. Thus, indeed, all square brackets are 
symmetric. Moreover, for B = A, we see that [A, A]i2 > and that the zero is only reached 
if A is conserved. 

Combining the results, we find that the brace product ()E1|) is symmetric. Moreover, 
we find that the brace product of quantity A with itself is nonnegative, {A, A} > 0, and 



vanishes only when A is conserved. As the brace product has features o: 
product, a version of the Cauchy-Schwarz-Buniakowsky (CSB) inequality 



a_ pseudo-scalar 
holds, 



{A,A}{B,B}>i{A,B}f . (E4) 

All the results from this Appendix remain valid, in an obvious manner, when the brace 
product ()E1|) is generalized to the pairs of tensors of the same rank associated with the 
particles, when requiring that the tensor indices are convoluted between the two tensors in 
the brace, as e.g. in (jHUj) . The positive definite nature of the brace product is important 
in ensuring that expressions for transport coefficients, obtained in the paper, yield positive 
values for the coefficients that in this case represent a stable system. 



[1] J. M. Lattimer and M. Prakasli, Astropfiys. J. 550, 426 (2001). 
[2] F. Rami et al., Pfiys. Rev. Lett. 84, 1120 (2000). 



35 

[3] G. F. Bertsch, S. Das Gupta, Phys. Rep. 160, 189 (1988). 
[4] J. Aichelin, Phys. Rep. 202, 233 (1991). 
[5] S. Tomonaga, Z. Phys. 110, 573 (1938). 

[6] V. M. Gahtskii ct al., Yad. Fiz. 30, 778 (1979) [Sov. J. Nucl. Phys. 30, 401 (1979).] 
[7] P. Danielewicz, Phys. Lett. 146 B, 168 (1984). 

[8] R. Hakim, L. Mornas, P. Peter, and H. D. Sivak, Phys. Rev. D 46, 4603 (1992). 
[9] R. Hakim and L. Mornas, Phys. Rev. C 47, 2846 (1993). 
[10] P. Danielewicz, Acta Phys. Polonica B 33, 45 (2002). 

[11] S. Chapman and T. G. Cowhng, The Mathematical Theory of Non- Uniform Gases (Cambridge, 

New York, 1964). 
[12] E. A. Uehhng and G. E. Uhlenbeck, Phys. Rev. 43, 552 (1933). 
[13] E. A. Uehhng, Phys. Rev. 46, 917 (1934). 
[14] E. J. Hehund and E. A. Uehhng, Phys. Rev. 56, 818 (1939). 

[15] R. L. Liboff, Kinetic Theory, Classical, Quantum, and Relativistic Descriptions (Prentice Hah, 

Englewood Chffs, New Jersey, 1990). 
[16] L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Vol. 6 of Course of Theoretical Physics 

(Addison- Wesley, Reading, 1959). 
[17] S. R. De Groot and P. Mazur, Non- Equilibrium Thermodynamics (North-Holland, Amsterdam, 

1962). 

[18] G. Baym and C. Pethick, Landau Fermi-Liquid Theory (John Wiley k, Sons, New York, 1991). 

[19] V. Greco, M. Colonna, M. Di Toro and F. Matera, Phys. Rev. C67, 015203(2003). 

[20] S. J. Yennello et al., Phys. Lett. B 321, 15 (1994); H. Johnston et al., Phys. Lett. B 371, 186 

(1996); H. Johnston ct al., Phys. Rev. C 56, 1972 (1997). 
[21] L S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 
New York, 1979). 



36 



TABLE I: Diffusion coefficient Dj obtained witfiin different orders of calculation, using experi- 
mental np cross sections, at sample densities n and temperatures T in symmetric nuclear matter, 
for species-independent mean fields. The numerical errors of the results on Dj are indicated in 
parenthesis for the least-significant digits. The last two columns, separated by the 'it' sign, give, 
respectively, the relative change in the result for the highest calculated order compared to the first 
order and the error for that change. 



n 


T 


Di 


Relative 






1^* order 2'^'^ order 


3^*^ order Change 


fm-3 


MeV 


fmc 


% 


0.016 


10 


0.29949(15) 0.3055(12) 


2.0 ± 0.4 


0.016 


60 


2.3891(18) 2.390(14) 


0.0 ± 0.6 


0.16 


10 


0.27964(21) 0.2800(29) 0.2809(25) 0.5 ± 0.9 


0.16 


60 


0.29591(24) 0.2965(19) 


0.2 ± 0.7 


0.32 


10 


0.4446(15) 0.4465(26) 


0.4 ± 0.7 


0.32 


60 


0.18187(15) 0.1827(13) 


0.5 ± 0.7 
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T (MeV) 

FIG. 1: Isospin diffusion coefficient Dj in symmetric matter, for Ui = 0, at different indicated 
densities, as a function of the temperature T. In the high-temperature limit, the diffusion coefficient 
exhibits the behavior Dj oc VT/n. Correspondingly, at high temperatures in the figure, the largest 
coefficient values are obtained for the lowest densities and the lowest coefficient values are obtained 
for the highest densities. In the low-temperature limit, the diffusion coefficient exhibits the behavior 
Dj oc n^/^/T^ and the order of the results in density reverses. 
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FIG. 2: Isospin diffusion coefficient Dj at normal density n = no = 0.16 fm^^ and different 
indicated asymmetries 5, for Ui = 0, as a function of the temperature T. An increase in tlie 
asymmetry generally causes a decrease in the coefficient, as discussed in the text. 
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FIG. 3: Thermal conductivity k in symmetric nuclear matter, at different indicated densities in 
units of no, as a function of temperature T. The conductivity increases as density increases. 
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FIG. 4: Shear viscosity t] in symmetric nuclear matter, at different indicated densities in units of 
no, as a function of temperature T. The viscosity increases as density increases. 
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FIG. 5: Mean-field enhancement factor of the diffusion coefficient in symmetric nuclear matter, 
R = Di{Ui)/ Dj(Ui = 0), at a fixed density n, as a function of temperature T. The solid and dashed 
lines, respectively, represent the factors for the assumed linear and quadratic dependence of the 
interaction symmetry energy on density. The lines from top to bottom are for densities n = 2no, 
no, 0.5 no and 0.1 no, respectively. At normal density the results for the two dependencies coincide. 



